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ABSimcr 



Whenever active or passive devices are ijrcliided: iir air 
electronic circuit that is to be analy^red by a crcnrputex., 
appropriate models for these devices imiSEt: bee deevecLjocpesd:, . 

The majority of the models that have beear. adapted: fixrr anarr- 
lyzing a circuit on a computer have negrLected: the: efdeectt off 
ambient temperature on the element* Tenrpeera?ture varraflen- 
has been shown to have hazardous effects on. eLectroiilic 
components . This paper deals with the mu dn-fn-cat n -mr of s.omee 
existing models to include temperature effects and: showss thee 
results of such modifications. It is done for. six, basic: 
electronic components (i.e., resistor^ capacitor, inductor, 
diode, bipolar transistor and field-effect* trans is tor) . 
Suggestions for extension of this work are also included. 
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I . imiEGrDTj’crTON 



A. GENERAL 

It is necessary far tke electrajiic: diesslgpT. eaigdireerr to: 
predict and identify tE.e reaFpcms^e of edleacltnDcrrix: sysfteorss to: 
temperature effects caused, by tftie- djeEv.iinse iitsjedjE orr tttee 
environment. Therefore ^ in. tiLe study- erf' semijCDiTdiictorr 
electronics the desirability a£ cunrpxLter sdjniilatlnn' become s- 
readily apparent. The equations desca±b±ng: sejiriJcnnductor: 
behavior in all states are compLex: and: nocniineaT:; : hence^ . 
solution by manual means is dafTieuXt ev.ecn. when', many • 
simplifying assumptions are made.. 

B. BACKGROUND 

1 . Circuit Analysis by Computer 

When a computer is used to analyze an electronic 
circuit, a sufficient amount of information concerning the 
circuit parameters must be supplied and device terminal, 
relationships in the network equations describing the 
circuit must be determined. This is not a difficult' task 
for passive components such as resistors and capacitors. 

Active devices may be described to computer analysis pro- 
grams in terms of models which, mathematically simulate 
the behavior of the devices.. It has become increasingly 
important to develop accuirate models of these: devices: for' 
computer programs used in design and analysis of ' such circuits. 
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Analyzing the netwarfes requires Idratt tdte cieaiia/'.edd 



equations must be in some paEirticuJiaiir jEcnnir. srd: tiiiedirr parnmee- 
ters specified properly by tbe user. Thre? ciircuiit: edj einen t35 
used in a network could include lumped resisdrances induc- 
tances, and capacitances, cumrent: and •vrolLtageE sa3nrrcES5 
(dependent and independent) , and acrCi'vre d^ev.iice? iiimiedLss (iL,ei.., . 
diodes and transistors). 

In the application of the Computer Aided Ana?iyS'is= 
Programs the need often arises for a circuit elemenr. that: 
is not included in the list of standard elements^.. Thee 
necessity of using linear but nonstandard', eiementss (such: 
as small -signal transistors or ideal transformers) arises>' 
in many programs. In the case of linear nonstandard ele-- 
ments , general equivalent circuits have often been developed 
by others; they can be applied to programs. 

2 . Simulation of Circuits with Active Elements 

Active networks contain dependent current and 
voltage sources or negative resistors.. When we attempt to 
analyze an active circuit by computer we encounter some 
problems. The first one is to find a set of mathematical 
equations to relate the current and voltage characteristics 
of each device. These equations can be obtained by one of 
two ways ; they can be derived from mathematical description 
of the physical phenomena occurring during the' operation 
of the device or might be obtaiued by erapiricall ohservations.- 
of the behavior of actual device. The degree of the 



8 



prediction, over a wide ircErng:® gd£ GqpTeai^ttdioTn aimiitt axttrall 
behavior of the device, detemniiiTesE thee varliiee ocff tliesBe 
equations . Secondly the panranretenrs^ irr. thessec eqnatlbnss must 
be known or measurable or the required irrforjrratioir. should, 
be derivable from temiiaaL meaESureraeirtsf.. Tfiee octherr fartorr 
is the tolerance of the paranretersE .. 
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II. DEVICE -MOD EEL INC METHODS' 



A. INTRODUCTION 

The objective of madelLijig: is ter prroiCEsesi; feoiir thisr 
fundamentals, that is, frcnir. matCTXsl pxoqienrtriess andfirror?- 
phology of a distributed device te s nrerdi^ ass po:eseiTted_ 
below : 

(1) To obtain systematically a lumped,. netwoTi—like: 
model for a distributed device,, a) emre shmilji wrilBr 
equations which are distributed airaia^s toe oerree oif 
Kirchhoff’s two network laws, the currenti-nade: lew 
or the voltage- loop law. Since the desiired: modell 
is network- like , it must obey the two laj^s- of-' 

Kirchhoff. b) One should remove, by some method 

of approximation, all aspects of the mathematical 
description which can be traced to the distributed 
nature of the device. 

(2) To obtain models which describe the transient 
performance of a device, one might take advantage of, 
if possible, knowledge of its steady-s-tate properties. 
The reason for this is, of course, that many aspects 
of steady- state behavior are much easier to deter- 
mine than their transient counterparts. 

In practice, a model is rarely perfect... Thus = we. 
should search for a model whose elements provide the best* 
combination of linearity and frequency independence. 
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In this section the £oer foasic a-:ppiroaKiTesE wiiUll bbee 



presented for developing equivalent ciroidits forr edi emeaits f 
that are not standard in computer programs. Theset approach es.- 
can be used individually or can be combined to form more 
complicated device model. 

B. PIECEWISE LINEAR MODELLING 

In piecewise linear modelling technique neTrliTreair 
elements can be analyzed by replacing all the nonXinear 
components with equivalent circuits ^ constructed wiith. 
the linear elements recognized by the computerr progranr.. 

The nonlinear characteris tics are averaged oveir the range? 
of interest and represented approximate ly^ by lineariized 
characteristics. For devices passing from one operating 

4 

region to another (such as a transistor moving from the 
active region to saturation), a different, linearized 
equivalent circuit can be proposed for each region. 

C. MATHEMATICAL MODELLING 

Those circuit elements that are not standard (defined 
as a circuit element that is not included in the list of 
built-in models of any computer library) can be modelled 
using only the mathematical relationship which describes 
their terminal characteristics. The procedure to follow 
in developing equivalent circuits using mathematical model- 
ling technique is quite straightforward.. En. many casus-,, 
the modelling problem can be reduced to the synthesis of' 
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a simple RLC network scinie times reqxidiiriin:g: tdree usase of: 
dependent current sources. Im otlteir ccssess,, ^ lliittttLiee 
inventiveness is useful. 

The equations describing the conrponents musrt be:- Unear 
since only linear elements are used in the- dev.eLapineiit: of; 
the model. The availability of L,. and C edementss inr. 
the construction of the equivalent circiiiit: permits? the 
modelling of linear differential equa-t±ans wiith: consitant;: 
coefficients . 

D. PHYSICAL ANALOGY: PEECEW:rSE LINEAR MODELLING 

The modelling of a nanlinear elenienl: wiith: arr. ali^biraii:: 
voltage-current relationship (such as ncmlinear- resdsitance:) 
was discussed in Section B. In th±s section the techniques 
for modelling a nonlinear reactive element with a non- 
algebraic voltage-current relationship is discussed (such 
as inductance or capacitance). For better understanding 
one may refer to examples given in Ref. 17. Since non- 
linear inductance and capacitance are algebraic functions 
of magnetic flux and charge, respectively, which are not 
recognized by the program, a reasonable solution to this 
modelling problem is to find a physical analogy which will 
allow flux or charge to be treated in those units recog- 
nized by the program. 

The physical analogy acts as a black' trox which converts 
the physical relationship into terms which are acceptable. 
to the program as shown in following algorithm. Therefore, 
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elements such as nonlinear indntcltaiices airdi csi^jmiiitajTCESs canr. 
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Figure 2-1. Physical Analagy Concept: 
in Physical Conrporteiit ModelLing;.. 



be modelled if an analogy betwreen the chnractreniisLtd'.cs? to: ber 
modelled and a voltage-current relationship acceptnbLe' to 
program can be established. 

E. PHYSICAL ANALOGY: MATHHLATICAh MODEhLXNC 

Almost all of the Computer Aided Circuit: AnaCysis: 
Programs perform calculations in terras of voItagcS' and 
currents only. The parameters, such as chp^rge and magnetic 
flux, have not been included in the programs and cannot be 
used directly. If such conditions exist, it is necessary 
to develop methods for expressing the new parameters (charge, 
magnetic flux, etc.) as functions of voltage and current or 
to redefine some new elements in addition to the elements 
included in the standard list that will perform the required 
analysis in terms of voltage and current. 

The solution of mechanical and thermal problems by 
electrical analogy techniques is an example of this- kind. 
Establishing the analogy is the first step in the dexiva- 
tion of an equivalent circuit. The second strep rS' linking- 
the derived analog network to the remainder of‘ the: model' 
which includes the element terminals. The element terminal 
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characteristic must always be defimeci in ttesmirs crE uodltaaggec- 
and current. The charge-control transistrcnr mmieL iiSE a 
good example of physical analogy techniques in madelllliiTg: 
the other charge equations which govern the transistror 
response. See Ref. 17, Appendix A. 2. 
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III. TEMPERATURE EFFECTS 



A. TEMPERATURE EFFECTS ON RESISTOR 

The resistivity p of a conductor is the ratio of the 
electric intensity E to the current per unit cross-sectional 
area. 

p = E/(I/A) 3.1 

The resistivity of all conducting material is affected by 
its temperature. A plot of resistivity versus temperature, 
for a metallic conductor, is given in Fig. 3.1. The curve 
may be satisfactorily represented by an equation of the form 

P = Pq ^ glT * § 2 ^^ * * 3-2 

where p^ is the resistivity at 0° C, g^, g 2 , ..., are 
constants depend on the characteristic of a particular material. 
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For temperatures which; are rccrt: tcrcr. g:ceaitl t^ies ttaainiss iir. 





3c. 3? 



The temperature coefficient definitiaEL a£ resis^t±v.iity iis- 



Since the resistance of a giren cCTEDciu;c:t:cn:- is 

to its resistivity the following equation can. tre wriitt^n 



where Rq is the resistance at 0® C and R is thee neshisstance; 
at T® C. 

The resistivities of nonmetals decrease with increasing 
temperature, and their a*s are negative.. Electrolytes' also. 



have negative a's. 

B. TEMPERATURE EFFECTS ON INDUCTOR 

Temperature variation is usually accompanied by a change 
in inductance, resistance and self-capacitance of a coil. 

The last two types of change are not important in compari- 
son with the change in inductance. The formula for single- 
layer solenoid given by Sturley [Ref. 6] is 



a '= gj/pQ = Cp-pgl/ChgT) 



Si.4: 




R = R^(l + aT) 






L = (r^N^)/(9r + 10£) 



3.6 



where I = length of winding 
r = radius of winding 



N = total turns in the coil.. 
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The change of the iintfaiC-tiainiEE fcir am iircxreeasjee iir, teurpearattiree 
of T® C can be stLowit ta he [[Hef.. 66]’ 

iLL = hL(A:T3 3:7’ 

where AL= change in inductance 

b = coefficient of iimeaEir ecspamsirmr. ocff thee conductor: 
AT = change in tempeir^ture.. 

Equation 3.7 is \raXid when, the cu3ecfficl:errts; off mdiail and: 
axial expans ion y and h^.,. ane app^ccxinraiteiy equal!. . 

It is possible to make the inductance: L_ independent: of : 
temperature by satisfying Fq:.. 36.. S,. the ratio: off axiaif to: 
radial expansion. 

hjh = Z + (9r/10£) 3f8' 

Xr IT 

By choosing the coil -former- and the conductor material 

4 

properly, the temperature coefficient for inductors can be 
minimized. 

C. TEMPERATURE EFFECTS ON CAPACITOR 

Temperature variation has two effects on capacitors, 
namely to change the dielectric permittivity and the dispo- 
sition of the plates. Some typical capacitance- temperature 
coefficients will be given in Chapter IV. Permittivity of 
the insulating material is considerably changed with tempera 
ture. Capacitances for parallel plate and cylindrical types 
are given with the following equations, 

parallel plate:: C = (A..e:)/d! 33 9:- 
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I 
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cylindrical: C = 



35.1DC' 



where area o£ plate (A) , and leng:t&: crE cmpncditmr ( jtO?,, aane: 
temperature dependent also, but this dependence can be 
negligible when compared with the deprendeuTce or£' p-ermitti-vity . 
In both cases capacitances are the functxon. of pTeriTii±h±iAity 
e. In addition to temperature dependence a£ rt tiss 
dependent of frequency, but the latter efifect iiS' a-- diiffecnent: 
story. 

Improvement of the manufacturing: techniques? mlhimlzjesr 
the temperature coefficient for capacitorrsF.. 

D. TEMPERATURE EFFECTS ON DIODE 
1 . Introduction 

There are two kinds of charge -carrying' particles'; 
one is negative (the free electron) of mobility and 

the other is positive (the hole) of mobility These 

particles move in opposite directions in an electric field 

E, but since they are of opposite sign,, tlie currents are in 
the same direction. Hence the current density J is given by 

J = (npj^ + pPp) e E 

where n = magnitude of f ree-electron concentration 
p = magnitude of hole concentration 
e = electron charge in coulombs. 

For an intrinsic semiconductor n = p = n^,. whera n^ is- the 
intrinsic carrier concentration.. The intrinsic carriex: 
concentration n^^ is temperature dependent, such as 
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n? = AT^exp(-E^Q/kT) . 

It has been shown experimentally that depends upon 
temperature also as pointed out below 

° ^GO • 

where E^q is the energy gap of the semiconductor at 0° K 
and q is a constant (q = 2.4 x 10 ^ eV/°K for silicon and 
q = 3.9 X 10 ^ eV/°K for germanium).* Such a large change 
in conductivity with temperature places a limitation upon 
the use of semiconductor devices in some circuits. 

2, Temperature Dependence of Parameters 

For a p-n junction the current I is related to the 
voltage V by the equation 

4 

I = I^[exp(V/mV^)-l] 3.11 




*Strictly speaking, q_is not a constant. For example, for 
silicon q = 2.4 X 10 eV/°K only when T ^ 200° K. For 
germanium, q = 3.9 x lO""^ eV/°K only when~T ^ 150° K. 
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J 

I 



is called the reverse satunratiiaii curirejiUi amK gijv.egj- }±y 
the equation 



= K t" expC-V(^/mV^) 



r.UUa 



~ ff-7~SSr vocLtS' 



for Germanium: m = 1 n = 2 

for Silicon: 

From experimental data we find trhnt the rev^s-e ssatiiratlpn- 



m=2 n = l.S vaJlts 



current increases approximately seven per cent: perr C.. 
Since (1.07)^^ = 2.0 reverse saturation, current: approximate: 
ly doubles for every 10° C rise in tempemture.. i^ the: 

"volt equivalent of temperature" defined; by 



= T/1I600 . 

Another important parameter for diode is diffusion, or storage: 
capacitance and given by * 

Cjj = T I/m i.12. 



where t is the diode time constant equals the mean lifetime 
of minority carriers. 

Reverse leakage resistance also has some temperature depen- 
dence; since it can change depending on the manufacturing pro- 
cess, no generalized temperature coefficient can be siated. 

It should be determined for each family of the diode. 

The temperature coefficient of the diode paramete.rs can 
be derived in differential form from Eqs. S. 11 and 3.11. a 
by taking the derivatives and making- some approximations 
to yield the following equations: 
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35 13: 



dV/dX = [V-(^Qo nnrV;j.)I/T 35114- 

at fixed current. 

E. TEMPERATURE EFFECTS UE TEANSE^mR 
1, Introduction 

A junction tranststax- ainatfrtSr ocE a: sdUicDn' orr 
germanium crystal, in w&icE a Layex of n.-type: orr p7-type: 
material sandwiched between two layers^ ocE p—type: orr n'-r type 
material, is called a p-n-p- cnr n-p'^m trransxsstocrr. , Forr ar 
qualitative understanding of the fcmiL of' the input: and: out 
put characteristics, one can consider that the. transistor 
consists of two diodes placed in series "hack to back," I 
has four operating regions defined as follows : 

Region (a) : Cut-off region (emitter reverse biased, 

collector reverse biased) , 

Region (b) : Normal transistor operation (emitter 

forward biased, collector reverse biased) , 

Region (c) : Saturation region (emitter forward biased, 

collector forward biased) , 

Region (d) : Inverse transistor operation (emitter 

reverse biased, collector forward biased). The 
input- output characteristics of the transistox in 
different modes are given in many texts.. 
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2. Temperature Dependeiice erf Fajimirgters- 



All of the parameters of orterresst: itr. computerr 
application concerning transistors chang;e wi±h". tenTperature. 
Some of them are very sensitive to teanpeLrature.- changes^ 
while others may be assumed essertiiaLLy craisstaaitt inn thee; 
analysis. The most temperatuxe sensciit±v.ee parameters 5 are e 
saturation currents, and (^^rrentn g^irr. factersf, , 

and Oj, and leakage resistances,. and: R^,. Prarameters; 
like time constants, x^, capaciitancns- C^,, C^. pTopor>- 

tionalities, m^, m^, change only slrgh-tlly' overr thee 
temperature region of interest as high: ass L25?' CZ. 

(a) Reverse saturation currents : The emitterr 

junction diode saturation current and coTiecterr junction" 

diode saturation current are the most temperature sensi- 

tive parameters. Suppose that the reverse saturation 
current at room temperature (300° K) is given by Eq. 3. 11. a; 
it can be written at temperature T + AT, as follows: 

I CO = h expC-Eco/kT) 



IcoT ” exp[-Eco/k(T + AT)] 



3:. 15 



equivalently 

= K^CT + AT)" exp[-Eco/k(T * ^T) ] 3.16 

Assuming AT << T and expanding the Eq .. 3:. .IIS' in Taylor series, 
the follov\ring can be obtained: 



Z2 
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^coT = mmi = n^[[]i Hm"j\] s^iiz' 



After numerical substitution it tunis out: to: tre 



IcoT = 0.157 AT) 



where X = 0.157/®C and is defined as tJxe teaipm:atiir:ee^ 
coefficient for silicon transistor. Ectir grreatrer:- atrcuraxy , . 
over larger temperature ranges^ tlie equa>t±an:. iinToiiiding: 
the exponential factor directly should be used., 

(b) Current gain factors r The current: giaiTr. factors; 
in normal and inverse modes for common, base conffguratton: 
are given below: 

“n ' ^ ‘‘FE^ “l “ ^ ’"EEr^ 

where hp^ is normal mode, common emitter current: gain,, and 

4 

hpEi is inverse mode, common emitter current gain. and 

aj do not change appreciably xvith temperature but (l-Uj^) or 
(1-aj) plays an important role in the following formulas: 



FE 



= aj^/(l-a^^) and h 






FEI 






For example, the variation of common-emitter current gain 
hpp, with temperature normalized to unity at T - 25° C 
for the 2N524 transistor is shown in Fig. 3.5. 

(c) Reverse leakage resistance: Reverse leakage 

resistance is temperature dependent to some extent. No 
generalized temperature coefficient can be stated since it 
varies, depending upon the process of manufacturing.., Thene-- 
fore, it should be determined for each family of transistor'. 
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F. TEMPERATURE EFFECTS ON FIELD-EFFECT TRANSISTOR 
1 . Introduction 

Field-effect transistor (FET) has two primary 

operating regions, namely the triode region in which the 

device behaves as a voltage-var iable resistor, and the 

pinch-off region where the current is determined mainly by 

the gate voltage. The two regions are separated by the 

locus of pinch-off points, i.e., the points where = 

V - Vnc» Vr.c » and are defined drain to source 

p (j O L) O ir Cjij 

voltage, pinch-off voltage and gate to source voltage, 
respectively. 

The operation of FET depends upon the flow of majority 
carriers only. Therefore, it is a unipolar device. Since 
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the maj ority- carrier cminreiit: (fedreases w.iitiii tteampEaiadrLLne^ . 
the troublesome phencraieiiaii (TEe£^ 4^1 is mxtt ernimirrtmxdf 
with field-effect transistors. Tli±s pxopearty pmocv.iidjesf 
thermal stability for FET. 

2. Temperature Dependence o£ Paxaiirgtigrs 

Drain saturation current,. wltiicir. depeards? upon'' 

channel geometry, impurity density and: charger cnrrimr 
mobility, is defined to be the channel current,. : 
at Vgg = -Vp. Although is not greatly dependent on 

Vpg, it depends on temperature heariily.. leraperature: can*', 
influence the drain current in two ways.. HiirsEtly-,, as changee 
in mobility causes a change in resistance cr£ channel.. 

-Secondly, a change in temperature causes a. change: in. siirfhce^ 
potential. The change in sur£ace potential works- in' opposi-- 
tion to the change in mobility. Thus, change in surface 
potential constitutes a built-in gate-voltage compensation 
of the mobility change in the semiconductor. The effect is 
dependent on the D.C. gate bias so that at a particular bias 
the compensation may be optimized. 

Another quantity of interest is the transconductance. 
The transconductance, g^, has the same temperature variation as 
Ipg does. The principal reason for the negative temperature 
coefficient of is that the mobility decreases with 

increasing temperature. 

As pointed out before, the major el'fects o£,' tempera- - 
ture on a FET or alteration o£ the average channel conductivity 
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Figure 3-4. Normalized Ij^gg versus Tempera- 
ture, at 25® C 



and the gate-channel junction barrier potential. Since 
conductivity decreases with increasing temperature, the 
drain current decreases as the temperature rises. On the 
other hand, barrier potential decreases with increasing 
temperature, causing the magnitude of the total gate-source 
potential to be reduced, hence increasing the drain current. 
By suitably choosing the operating point, it is possible to 
make these effects canceled, producing what is essentially 
a zero temperature coefficient operating point. 
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IV. COMPUTER MODELS 



A. INTRODUCTION 

Almost all the computer analysis programs are designed: 
to solve the linear and nonlinear simiiltan.eaus equ;at±OTts 
which characterize the mathematical models us^ed tn: presdxct: 
circuit response for an electrical/temperature eTTviranment:.. 

The user codes the problem by entering the circuit topalogy 
and parameters. Programs contain a standard (buiiit-in.) 
component model library. Non-standard models may he- cnn- 
structed by using standard models. FORTRAN algebraic 
equations and control statements may be coded to modify 
standard model parameters, augment the standard models, 
create nonstandard models, create additional time functions, 
temperature functions, etc. The circuit is entered into 
the computer program by using the program standard part 
models and specifying the circuit topology. 

In this paper by '’conventional (or non-standard) model," 
it means the model of a device which is isolated from all 
environmental effects, and by "modified model," it means a 
model with temperature effects taken into account. In 
addition to this task, some additions or some simplifications 
will be done when necessary, depending upon device characteris 
tics . 

In order that electronic components may be handXed 
quantitatively, a set of mathematical relationships must, be 
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derived relating their terminal characteristics. Deriva- 
tions will not be attempted in this paper since they are 
documented elsewhere [Refs. 2, 3, 4, 7, 18, 22]; only the 
results will be presented. 

B. RESISTOR MODEL 

1 . Conventional Model 

The resistor model is defined by the following 
current equation 

iR = V^/R 

where Vj^ is the instantaneous voltage across the resistor R. 

^R 

> ^ 

T, T. 

R ^ 

4 

Figure 4-1. Resistor Element 

Since the element is bilateral, the choice of which is con- 
sidered to be terminal 1 and which is terminal 2 depends on 
the direction of current flow (from T^^ to T 2 ) . This element 
has just one parameter: the resistance value R. 

2. Modified Model 

Temperature dependence of resistor was discussed in 
Chapter III, but the Eq, 3.5 can be rewritten in slightly 
different form such as 

R^ = R ± AR 4.1 

where AR is the amount of resistance change due to tempera- 
ture change and equivalent to 
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AR = ±aR(AT) 



4.2 



The sign depends on material properties. Assuming resistor 
is linear, the following modified resistor model for 
computer can be used. The determination of AR requires the 
parameter a, which can be found in many manuals or reference 
books . 



Ti* 

Figure 



^ -"R 



R AR 

4-2. Modified Resistor Model 



T 



2 



C. CAPACITOR MODEL 

1 . Conventional Model 

A capacitor, shown in Fig. 4.3, is a 'bilateral 
element with terminals T^ and T 2 , having a capacitance C. 
The current flowing through a capacitor is given by 

i^ = C[d(V^ - V2)/dt] 



Ti* 1 ^ ^ ^2 

Figure 4-3. Capacitor Element 
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2. Modified Model 



The more realistic capacitor has a shunt leakage 
resistance and a series resistance as shown below that are 
identified with and R 2 , respectively. 




Figure 4-4. Capacitor Model with Resistive Components 



Temperature dependence of capacitor was discussed in 

Chapter III. Very small series resistance ,(in thfe order of 

0.1 ohms) and very large shunt resistance (in the order of 
9 

10 ohms) can be assumed constant with temperature change. 
Some typical capacitance- temperature coefficients are given 
in the following table and many others can be found in the 
references given in the footnote.* 



Table 4-1. Cap 
Insulating Material 
Varnished Cambric 
Synthetic Resin 
Enamel 
Ceramic 



acitance- Temperature Coefficients 

Capacitance- Temperature Coefficients 
+2100 parts in 10® per degree Centigrade 
+1600 parts in 10® per degree Centigrade 
+470 parts in 10® per degree Centigrade 
+100 parts in 10® per degree Centigrade 



*(1) Handbook of Chemistry and Physics; (2) International 
Critical Tables Vol. 2, McGraw-Hill; (3) Tables of Dielectric 
Materials, Vols. 1-IV, MIT, January, 1955. 
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Defining a change in capacitance as AC the modified capacitor 
model can be obtained shown in Fig. 4-5, in which AC should 
be calculated and added to the model with the aid of tables 
(given in the footnote on page 30) and the Eqs . 3.9 and 3.10. 
The change in capacitance is very important if the capacitor 




Figure 4-5. Capacitor Model with Temperature Effect 



is an element of resonance circuit because a small change in 

4 

capacitance may cause a large drift in resonance frequency. 

D. INDUCTOR MODEL 

1 . Conventional Model 

The relationship for current is given by 

i^ = (1/L)/ (V^ - vp dt 

where L is the value of inductance in henries and the dif- 
ference, (V^ - V 2 ) , is the voltage across the element. The 
inductor is diagrammed in Fig. 4-6. Since the inductor is 
bilateral, the choice of which to be considered T^^ is com- 
pletely arbitrary. This element has two parameters: induc- 

tance L, and Q (of the coil). Many applications of inductors 
have practical equivalent circuit that also includes the 
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Figure 4-6. Inductor Element 



quantity Q--not just pure inductance. Frequency dependent 
r is given by 

r = Q.ZirfL 4.3 

where £ is the frequency at which the analysis is being 
performed. The significance of Eq . 4.3 is best seen in a 
sketch of resistance versus frequency. The range for the 
Q model is that r becomes unrealistically small as frequency 




becomes very low, and at high frequency the model would show 
r approaching infinity which is also unrealistic. 
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2. Modified Model 



Temperature dependence of inductance was given in 
Chapter III with Eq. 3.7. Additionally D.C. series resis- 
tance and the shunt leakage resistance should be involved 
into the model, which are assumed to be independent of 
temperature. The schematic which takes all the effects 
explained above into consideration is shown in Fig. 4.8. 



AL 




To determine AL one should first determine the parameter 
b in Eq. 3.7. The value of b can be found in the handbooks 
specified in the footnote in previous section. For better 
understanding of importance AL, look at the example given 
below. 

Example: For a tank circuit, f^ = 1/ (2 tt/LC) . If Af is 

the change in oscillation frequency due to AL change in L, 
the folloxving relationship can be obtained after proper sub- 
stitution and expansion by the Binomial Theorem where only 
the linear term is retained: 



Af 

T “ 




L^ 
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The value of b far Cdrjrpeir is II..6 d x>. TOC ^ paa:?t:ss ppar ®’ CC. Thus 5 
the ratio change in firequHixc^- (fuje; tCD thee urrhliidhanedc e^xpansionr. 
of a copper coil is 8 x 10 paeirts peEir ®'' C,. ocrr furr AT7 == 3D.® CC 
the ratio change of frec^uency is 24-0 parts in T0°.. Firequency 
variations due to ceil expsErrsian! oh 24TC ^rdf 2I4-OT. hzz att oscilr - 
1 at ion frequencies a£ I and TO' Mhr,. rr^qp^c±±v.eiy-, . woulddbee 
expected under these ccmfitions.. 

E. DIODE MODEL 
1. General 

Any diode can be nrcnEeTTed: using: thee ffxlluwiirg: approach . . 
The quantities of interest are eTectricaell chaTnctexisti cse 
which can be observed at the externaT terminarse (ii.e^., .aa 
”black-box*' approach) . This type of nradeO'ing. requires.' 

4 

finding a mathematical expression to represent each electri- 
cal characteristic. In this paper this method will be 
emphasized. Physical constants will be used whenever possible, 
but if a characteristic cannot be related physically, an 
empirical description will be used. 

Based on the operating range, diode models may be 
classified into two groups: '’small -signal” and "large-signal” 

models. The first type consists of finding an operating 
point and assuming that all responses are linear about that 
point. Therefore, this model is valid over only a limited 
operating range, namely, where the response characteristic 
is linear. Large-signal model is one which can- be. used„ over, 
a wide range of operation. By its nature, if is the best 
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suited model for use in geneiral ccrntpuiteir arr^^srsE.. Tlregn^ 
are three most widely knawn mxKiels: (T) Khrgrrs^-M fanlll'! ID^ 

(2) Beaufoy-Sparkes charge- contra 1 model ['He£.. Ill] and: 

(3) Linvill lumped model [Ref. 12]. In this paper Eirers- 
Moll model will be examined in nnore detail sni neaiESSsary 
modifications will be given. The reason forr thisE preiiecnencer 
is this model describes the "hlack-hox.^’’ hehaEv.iorr with', redas-- 
tively simple mathematical equations. 

2. Model without Temperature Effects 

The current equation for diode was g^ven, bnefbrec ihr. 

Eq. 3.11. Another diode current equation, w.illli he oirtaineii 
by discussing the various factors in the semiconductor: dinde: 
which have to do with its electrical performance.. 

With the application of an external reverse bias, 
the p-n junction width will expand because more donors and 
acceptors will be uncovered. This increase in uncovered 
charge with applied voltage may be considered a capacitive 
effect, named transition-region capacitance.. 

If bias in the forward direction, the potential 
barrier at the junction is lowered and holes from the p-side 
enter into n-side. Similarly, electrons from the n-side move 
into p-side. A change of voltage on the p-n junction will 
alter the amount of stored carriers. This, then looks like 
a capacitor to the accompanying circuit, named storage or 
diffusion capacitance. 

By taking the leakage resistance into account,, the follow 
ing diode current equation can be written as 
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I + E/Rj^ 



Cp(dE/dt) + C^(dE/dt) 



4.4 



Id = I ^ ^ - 

These previously derived properties are now all combined to 
provide a composite equivalent circuit, or computer model, 
of the p-n junction diode, shown in Fig. 4.9. 
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Figure 4-9. Diode Computer Model 



The definition and determination of the model parameters 
will be discussed in the following section. 

3 . Determination of Model Parameters 

Ideal model parameters can be read directly or 
extracted from the device specification sheet. Unfortunately, 
this cannot always be done. Most data sheets do not contain 
sufficient information to define the entire model. 
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In Fig. 4-9, Rg represents the s eraicaiiciuictGrrr feuUM: resEiiSc-* 
tance, the junction leakage resistance, the^ transciit±air 
capacitance, and I is an ideal current source whiich degcribres.- 
the forward current as a function of junction voLta:ge E.. It. 
is possible to write the Eq. 3.11 in slightly (Irffeergnt: focmr 
as given below: 



I = I^[exp(eE/mkT)-l] 4J..5c 

where e, k, and T are electronic charge, BaItzmamT.''g consetnnt:, 
and absolute junction temperature in °K, respectively.. I^ 
and m are constants which define the shape of th.e V-E dinrac:- 
teristics. For D.C. operation, the network equntiorrsF are::: 






Ij, = I * E/R^ 



4;..e 

4.7 



A plot of the static - 1^ characteristic for 
positive Ip is shown in Fig. 4-10. The current is plotted 
on a logarithmic scale to illustrate that over a wide range 
in the forward operating region, Vp = E. For determining 
the unknown model parameters I^, m, and Rg one method is 
to curve- fit data points on the forward diode characterrsiic 
to Eqs . 4.6 and 4.7, where Ip is assumed to equal I. Another 
approach to finding the parameters involves the same assump- 
tion but not the curve fitting. One may choose two paints 
on the linear portion of the curve in Fig. 4.10 and label 
them (I^, V^) and (^ 2 * ^2^* this range may be expressed 
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with following: 



" ^o exp(0Vj^/m) 



4.8 



where 6 = e/kT. 



Ijj (amps) 
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Taking the logarithm of both sides and substituting 



(I^, V^), (I 2 , V 2 ) , then solving for m gives 

m = e(V 2 - V^)/ln(l 2 /l 3 ^) . 4.9 



Substituting either data point and the value of m found in 
Eq. 4.9 into Eq. 4.8, the following can be obtained: 

= I^/exp(eV^/m) = l 2 /expCeV 2 /m) . 

The difference between the linear extrapolation and the 
measured point is the drop due to Rg. Solution of Eq. 4.6 
results in 



where E^ = (m/6) In [(I^ + I^)/I^] = (m/6) In 

4 




From 

ohmic 



the reverse characteris 



leakage resistance, 



tic shown in Fig 
can be found to 



4-11, 

be 



the 
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There remain two more parameters to be determined, namely, 
the junction transition capacitance C,j. and diffusion capaci- 
tance Cp in Fig. 4-9. A typical curve of junction capacitance 
versus reverse bias is shown in Fig. 4-12. 




The analytical expression that describes is 

= K/(Vj - Vp" ' 4.12 

where K, n, and Vj are junction capacitance proportionality 
constant, grading constant and junction contact potential, 
respectively (Vj is approximately 0.9 or 1.0 volts). For 
Vjj = 0.0, = K. Therefore, at any point (C^, V^) 

° 

Solving Eq. 4.13 for n results in 

n = ln(C^/C^)/ln(Vj - V^) . 4.14 . 
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Another way of expressing the Eq. 5.12 imay tee fibcillmiiTrg;: : 

= 6(1 + I }/in27rf 44. 5355 

u o 

where the parameter £ is diode cut-off frequency arrd: eA^a^Lua- 
tion of f can be done by Kuno's* meth-od.. Tke cuti-crEf' 
frequency is given by 

£ = 1/2ttt 

4. Modified Model 

The temperature dependence of diode was= exam±ne:d: iir 
Chapter III. Temperature coefficient of saiturntdioir cnrr.ejit: 
and temperature dependence of the diode •vroltaige werrer g^^en'. 
in Eq. 3.13 and 3.14. Equation 4.15 is express^ed irr. raone 
convenient form as 

Cjj = Glx/m = (ex/mk) (I/T) . ‘ (4. .16) 

For simplification purposes, making some definitions will be 
useful at this point such as follows: 

A = ex/mk 4.17 

B = eE/mk 4..I8 

C = Fractional increase of reverse 
current with temperature. 

Temperature dependent reverse saturation current increases 
exponentially as indicated by the following equation: 

exp[C(T - T^)] 4^..I9 

*H. J. Kuno, "Analysis and Characteriration of p-n- JUnctlbn 
Diode Switching," IEEE Trans, on Electron Devices,. Vol.. ED- 11, 
p. 8-14, Jan. 1964. 
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where is the reverse cnirrent at tearnpiearatiirg TT aaidd iiss 

the reverse current at temp era toare T^.. TUree Scodiitdricrr: cxff 
Eq. 4.19 for C turns out to be 

C = ln(I^^/I^3/CT - T^) .. 44.2D. 

For example 1N3605 silicon dioxie bas^ th^e £bdXawiirrg~ spffisrcEfirar - 
tions : 

I =0.05 mA at T = 25° cl 
o ° t 

j Applied Volta g:e = —3d volts. 

= 50 mA at T = 150° C J 

Therefore 

C = ln(50/0.05)/(150 - ZS) = 0:..ffEE/ °r .. 

Proceeding with concise form of Eq. 4 ..16,. the following 
relationship is obtained: •• 

Cp = A(I/T) . 4.21 

If we drop the unity in the parenthesis in Eq. 4.5 an approxi- 
mate current equation can be written as 

I = [exp (eE/mkT) ] = exp(B/T) 4.22 

where B is given in Eq. 4.18. The AT change in temperature 
causes the following modification for Eq. 4.22. 

If = IqT + AT) ] 

= exp(C.AT) exp[B/(T + A.T) ] 

= exp(C.AT) exp(B/T) exp[ -B (AT) /E(T. +- AT) ] • 

Definition of another factor such as 

D = exp(C.AT) exp[-B(AT)/T(T + AT)] 
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4.23 



provides a very short expression for I^, given below: 

= D.I . 

On the other hand AT change in temperature implies the fol- 
lowing modification on diffusion capacitance Cj^. 

= A.I^/(T + AT) = A.(D.I)/(T + AT) 

= (A.I/T) [T.D/(T + AT)]= . T.D/(T + AT) 

= F.Cj^ 4.24 

where F = (T.D)/(T + AT). Then adding and substructing the 
Cj^ to Eq. 4.24 yields 

CpT = + Cp - Cjj 

= Cp + Cp (F - 1) . 4.25 

' 4 

In similar fashion 



= D.I + I - I 



= I + I(D - 1) 



4.26 



The temperature treatment of more sensitive diode parameters 
was given by Eqs . 4.25 and 4.26. The changes due to tempera- 
ture for Rg and are assumed very small and neglected. 
Applying all necessary modifications one would get the com- 
puter diode model when temperature vatiations exist as shown 
in Fig. 4-13. 
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Figure 4-13. Diode Computer Model with Temperature Effects 



F. TRANSISTOR MODEL 
1 . General 

A transistor model in computer-aided analysis should 
meet the following properties: (a) The model must consist 

of circuit elements that are recognizable by the computer 
program being used. (b) The parameters used in the model 
should either be obtainable from published data or readily 
estimated. (c) The model should be valid to high frequencies 
as well as to low frequencies. (d) The model should function 
without regard to the circuit configuration (common base, 
common collector, or common emitter) in which the transistor 
is placed. 
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Note : A general model for a cfeTice does imtl; LeaGd; ttc 

the most efficient analysis since ceirCarm tairmfrlesaoDFiree eijeairerrts 
of the model may not actually be required by the pxcdrljenr ar: 
hand. For example, nonlinear or high frequency eleinentrs 
may not affect the circuit response gxeatly,, but: cmLid: 
easily affect the solution time significantly.. 

There are many models that hnve beear. pxaqnxsedi to: 
represent transistors at low frequencies.. Any of these; smalll- 
signal models is valid only when the transis-'tor operates* 
in the active region and produces neglig:ible disttrrtian. 
in the output signal. Fortunately the output chaructBr:is.t:acs5 
of transistors are reasonably linear over a fairly large" 
range, often allowing output swings of several volts with 
little distortion. High~f requency models :for the transistors 
are more complex than low^ frequency models, since additional 
reactive elements are included. 

2 . Model without Temperature Effect 

The transistor is a three- terminal device, fabricated 
by placing two p-n junctions in close proximity to each other. 
The current through one forward-biased junction strongly 
influences the current through the other. The normal mode 
of transistor operation occurs when the base-emitter junction 
is a forward biased and base-collector junctrion is reverse 
biased. A transistor may also be operated in the inverse- 
mode. In this case, the collector is forward* biased' and. the. 
emitter is reverse biased. Therefore,, transistor can be. 
considered as two junction diodes connected in series, of. 
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which one is forward and the erther is irevnsirseee fcrasagdi . TKec 
emitter junction and the collector junctxarr- currenlr ea:fLias- 
tions can be written as in single diode case.. Hy re^earrdjig: 
to Section E and taking the inverse mode current: gain. inte. 
account, the following current equations nra)^ lie ex:'prre£5ssaedi in. 
terms of diode parameters: 



^E ^EF ^^E^^LE^ C-pgCdVg/dt) +• Cj^g(d\^g^/dit) 

Ic = IcF ^ (V^/R^c) - C^c^dV^/dt) - S^(dV^/dt) - 4.2B 

where parameters and are normal and inventredl sliort:— 
circuit current gains between emitter and collectxirr,, capaxl:— 
tors C.pg and represent emitter-base dep letion-Layer’ and; 

emitter-base diffusion capacitances, respectively.. Whiie; Cij^- 
and represent the collector-base depletion and diffusion 

capacitances, and are the emitter-base and collector- 
base voltages, respectively. Emitter and collector reverse 
saturation currents may be involved into the junction current 
equations in such a manner given in the following equations: 





ijgF I F2 [ ® ^P C ® ) ~ 1] 


4.29 




I^P — I ^2 [ exp ( SV^/m^j - Ij .. 


4.. SO 


^ES 


and I^g are usually given in specification 


sheets , but 


the 


analytical expressions for I^g and are 


given as 




^ES " ^EO^^^ ‘ 


4:..3ZL 


and 


^CS " 


4 :. . 3 : 2 : 
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where is the saturation current of emitter Juimttiian'. wiith". 

zero collector current and is the satciraticuL current; ccf 

collector junction with zero emitter current. 

Now we are in a position to model the transistor 
device. The computer model uses two p-n junctians inxiderllLeEd. 
like the diode previously discussed plus the current g:^in. 
factors as shown in Fig. 4-14. The resistances R^. and: 
account for ohmic effects and can be approxinrated as constam 
in most instances. 

‘ The junction capacitance factors are: 



Emitter transition capacitance, 4^..3S 

and 

Collector transition capacitance, C,j.^ = (Vj - 4^,.J4- 

where n = 1/2 and n = 1/3 for an abrupt junction and a linear 
graded junction, respectively. C^q and C^q are determined in 
the same manner as explained in previous section (Section E) . 

Most transistor manuals give the normal and reverse 
common emitter current gains and The relatrion- 

ship between the four gain factors is given by 



^FE(N) 



"n ^FE(N)/^^ ^FE(N)^ 



TE(I) 



Uj/Cl - Uj) Qj ^FE(I)/^^ ^FE(I)^ * 



4.35 



The diffusion capacitances, C^p and are given by 



'DC 



^DE ®^EF *^e/^E 



and 



4:..36 



^DC ®^CF "^C^’^C 



4'..57 



47 



I 

I 



4 







where Xg and x^ are emitter and collector time constants, 
respectively. These time constants can be evaluated by the 
Kuno’s Method that has been mentioned in Section E. 




Figure 4-14. n-p-n Transistor Model for Computer Analysis 

Most of the transistor parameters are supplied by 
the manufacturer. However, if it becomes necessary to deter- 
mine the parameters. Refs. 18, 19, 20 and 21 are very useful 
for evaluation. 

3. Modified Model 

The temperature dependence of transistor parameters was 
discussed and some assumptions were stated in Chapter III. Addi- 
tional assumptions are given below: 

(a) and Uj are approximately the same at temperatures 

T and (T + AT) , 
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(b) AT << T, 

(c) Tg and are the same at tempreEiratuiiTes TT asudi (TT +- A!T) 

(d) Ief = Ies' exp (6Vg/m,g) I and the same apprroximatiixrr 
is for 

CF, 

(e) Temperature dependence e£ and: rss ass esqrLaaiieiii 

in Chapter III. 

Equation 3.17 was written far callectar jumrttixnr. ddbde:- 
saturation current, and a similar expressiioiT. canr. bres writ'tBn 
for emitter jucntion diode saturation current: as- 



EOT 



Ig^[l + X 






I 



)-I 



E0 



443SS 



where is emitter junction saturation current: at: 

temperature (T + AT) while I^q is at temperature T.. The: 
expression for ^^^s given by Eq. 4.31 for a specified 

temperature. If I. 
temperature (T + AT), the following expression can be written 



■cr^rr is defined as the value of I^rc 
EST Eb 



^EST ^EOT'^^^ ■ 



= [IeO ^ ^Eol/^1 - Vl^ 



^ ^ES ^E ^ES 



4.39 



where = A„(AT). Equation 4.L8 was the defining equation 
n h 

-for the factor B. A similar definition can be used: for 
emitter junctions as follows: 
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For the derivation of Igpj, first exprefss thee ttea::nr 
”exp[Bg/(T + AT)]” in a more conveixient fcTTi]ir. as fod-lxicwss: : 

exp[Bg/(T + AT)] = exp(B^/T) exp [-B^(AT) /r(T AT) ]’ 

= Lg expCBg/T) 

where = exp [ -B^ (AT)/T(T > AT)]. 

Thus for a change in temperature,. AT” „ thee fidJLowin.'g 
approximate expression for the forward emitter' cuirr.ent: is- 
obtained : 



hpT ■ ^EST 

= (IjS ♦ Eg IggJ [exp (Bg/T) ILg. 

= IggCl * Kg)Lg exp(Bg/n 

= Igp (1 + Kg)Lg 

4 

^EF 4.40 



where = Lg(l + K^) . 

Another definition given in Section E was Eq . 4.17 
which is repeated here for emitter junctions as follows: 

Af = eXg/Mgk 

Then can be expressed as = A^I^p/T. Therefore, 

emitter diffusion capacitance at temperature T + AT becomes 

^DET ^ ^E ^EFT^^^ 

= Np 4-. .4-1 



where = MpT/(T + AT). 
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As stated before, a IsEir^e (dxarrge? iiir. liecwearerr,, nray'-' 

reflect only a small change in a amf fm nrarry- irrsstaarness cair. 
be neglected in this treatment. Tlteirefam the: current', 
generator iri Fig. 4-14, where the temperature- varia- 
tion is not considered, can he written withr. 

temperature change (AT) , so that 

OjjIgFx " ^EF 44 4^- 

in which Eq. 4.40 has been applied. 

One may follow the same dertvationS' given, for the; 
collector junction and obtain the foUlow.ijrg: result3?:; 

Icft " ^CF 

C = N C 

^DCT ^DC 

ailfFT = ^CF * 

Same manipulation can be done for obtaining 
temperature-caused current generators and capacitors as 
in the diode case and then the following are obtained: 



^EFT 


^EF ^EF 


^EF 




^EF ^EF ^^^E 


- 1) 


^CFT 


^CF ^CF 


- 1) 


^DET 


^E ^DE ‘^DE 


r 

^DE 




^DE ^DE ^^E 


- 1) 


^BCT 


^DC Sc 


- I) 



Now the modified computer model for transistors can be drawn 
as shown in Fig. 4-15. 



51 



I 



I 






I 




Figure 4-15. n-p-n Temperature Effected T/ 

for Computer Analysis 

G. FIELD-EFFECT TRANSISTOR MODEL 
1 . General 

Although field-effect transistors 
used extensively at present, there is no r 
describe all such devices. One reason is 
processes causing various modes of operati 
another. There are two main types, juncti: 




' CT 



ll 
one 
; Lated 
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gate FET, sometimes called the "Metal-Oxide Semiconductor" 
(MOS) transistor [Ref. 13]. This would suggest the use of 
at least two different models. A piecewise linear model 
has been developed by Roberts and Harbourt [Ref. 9] that is 
designed for circuit analysis using the ECAP [Ref. 1] pro- 
gram. However, piecewise linear models have some deficiencies 
when applied to general circuit analysis programs because a 
piecewise linear model can be made to approximate the actual 
behavior quite accurately by using several regions. However, 
as more regions are added, more constants must be provided 
and the awkwardness involved in determining the correct region 
of operation. 

2 , Models without Temperature Effects 

It is possible to make analogy betwpen vacuum tube 
(such as a triode) and FET, provided that the cathode is 
replaced by source, the plate replaced by drain and grid 
replaced by gate. Neglecting the bulk resistances and taking 
the barrier capacitances into account, the following equiva- 
lent circuit can be obtained. Since the gate junction is 

Gate Drain 




Figure 4-16. Small-Signal JFET Equivalent Circuit 
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reverse biased, the gate- source resistance aarrdi thee 

gate-drain resistance are extremely lar^e,,. and: hearer 

have not been included in the circuit of Fig:^ 4^-L6c., Deeterr— 
mination of the equivalent circuit parameters will be- dis- 
cussed later. As pointed out before there aEree twee cutaeg^orlEse 
of field-effect transistors, JFET and MDSFET.. A^Sc thee names; 
imply, a 'JFET uses the characteristics a£ a reversEee-bl'asejd: 
junction to control the drain-source current,, whille in the. 
MOSFET the gate is a metal deposited on an oxide layer and 
is insulated from the source and drain. Both devices.- 
operate on the principle of a "channel'*' current controdled 
by an electronic field. The control mechanisms Ear the' two' 
are different, resulting in considerably different characteris. 
tics. The main difference between the two ^ is in the gate 
characteristics. The input of the JFET behaves like a 
reverse biased diode while the input of a MOSFET is similar 
to a small capacitor. Based on the above discussion, the 
models shown in Figs. 4-17 and 4-18 are obtained. 

In Fig. 4-17 the source end of the junction is less 
strongly reverse-biased or even slightly forward-biased, 
depending on the instantaneous value of It exhibits the 

voltage-variable capacitance and conductance properties of 
an ordinary p-n diode. A lumped model of this gate-to- 
source junction can be produced by using the standard built- 



in model for a junction diode to simulate these effects.. 
Transconductance and output conductance are simulated by the 
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Figure 4-17. Lumped Junction FET Computer Model 




Figure 4-18. Lumped MOSFET Computer Model 

nonlinear current source, which can be modelled by 

pointwise linearization as described later. 

As stated before, the gate of a MOSFET resembles 
a parallel-plate capacitor with a thin dielectric layer. 

The dielectric results in a very low gate conductance, and 
since the thickness is constant, the distributed capacitance 
is uniform along the channel and independent of This 

capacitance can be split between the source and drain ends 
of the channel and both lumped capacitors simulated by 
standard capacitor model. 
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For more general modelling (including large-signal 
behavior) by pointwise linearization, the functional relation- 
ship is not known in general but can be approximated emipri- 
cally for a given device by measurements, if the desired 
degree of accuracy so required. 




Transconductance g is the 

‘^m 



Figure 4-19. Pinch-off 
Region Transfer Curve 
Showing Pointwise 
Linearization 
Procedure 



GS 



slope of. the transfer 



characteristic curve. Fig. 4-19. 



At a point ^he 



tangent line is given by 



^DS ^m ^DS 

where 



4.43 



8m = - (^Gs/Vp)] • 

Ijjg in the pinch-off region is expressed as in Eq. 4.44. 



^DS ■ ^^GS^^P^^ 



n 



4.44 



where n is a constant for a given unit, whose value usually 
lies between 1.7 and 2.7 (for many units n = 2 is a good 
assumption), and is the drain current when v^g = 0. 

The determination of the necessary model parameters will be 
given in the following section. 
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3. Determination of Model Parameters 

Capacitance determination: For the model represented 

in Section 2, the common-source short-circuit input, output, 
and reverse capacitances are given in the following 
expressions, respectively: 







^ISS 


r + r 
^DG ^GS 










^OSS 


r + r 
^DG ^DS 






^RSS ^DG * 

Consequently, the model capacitances may 


be calculated in 


terms 


of 


the 


short-circuit capacitances as 










c 

^DG 


^RSS 










f 

^GS 


^ISS ■ ^RSS 


4 








c 

^DS 


Sss ■ ^RSS • 






If 


the 


capacitances 


^OSS» ^RSS 


are not supplied 


by 



the manufacturer the following measurement method can be used 
to determine the unknown values. Three terminal instruments 
are the most useful, since in addition to measuring the 
capacitance between two terminals, they are equipped to 
eliminate the effects of capacitance between these two termi- 
nals and a third terminal. This is shown in Fig. 4-20. In 
Fig. 4-20 the presence of and C 2 does not affect the 
measurement of C . In some instruments the "low" terminal 

X 

is grounded and the terminal connected to both and C 2 is 
called the "guard" terminal. Capacitance measurement jigs 
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are diagramed in Fig. 4-21. The limitations on this 
measurement are: (1) the gate signal should usually be 

limited to less than 200 mV, (2) frequency should be kept 
between 500 KHz and 5 MHz. Below 500 KHz the chokes must 
be large, and above 5 MHz series resonance and some other 
detrimental effects occur. 

Conductance Measurement: The pinch-off voltage, 

the drain saturation current, 1^33 » bhe drain conductance, 
gjj, and the transconductance g^ are four of the parameters 
required in the models presented in this section. For the 
pinch-off region, the transconductance is given by Eq. 4.45 



gni = (-2I^c;c;/Vp)[l - (V^^/Vp)] 



DSS' 



GS' *P- 



4.45 



if Vgg = 0 then g^ is maximum. 



or 



’max ^ ^DSS'^^P 



^P " ^DSS'^^max * 



The experimental circuit for determination of is given 

in Fig. 4-22. 



The determination of g^ can be done by a circuit 

as shown in Fig. 4-23, in which gj^ (micromho) = 100 V2(niV) 

when Vi = 1.0 volt. 

1 • * 



Measuring 

Ins trument 

Textronix 

Type 

130 L-C 

Meter 



^ — — 

High 

r -1 


L 


Low * 

Ground ^"^^1 


\ 

c 

X — 1 

] 



T 



Figure 4-20. Three Terminal Measurement 
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Figure 4-21. 



(a) 

(c) 



Cjss Measurement 
Crss Measurement 



(b) 

Jigs 



'Measurement , 



100 ohms 




Figure 4.22. Circuit for measurement of g^^^ 
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measuri 
can be 



Current Measurement: 



loss easiest determined parameter as far as the 

ng circuit is concerned. The following configuration 
set for determination of 




DSS 



10 



volt 
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4 . Modified Model 



As pointed out in Chapter tdie nr^jicnr eefffecirss of: 

temperature on a FET are to alter the aAreracgre cdrarmorL comiuc:— 
tivity, a, and the gate-channel junction barrxeTr pntentiall, 

Since a decreases with increasing temperatu,re,. the: effe ct: of:' 
this to decrease the drain current as the heniperrntiire r±S3es:.. 

On the other hand, ji|;| decreases wfth incneasEiing tenipeatatiLre:, . 
causing the magnitude of the total gate-scruarce pcrtentharll tn: be? 
reduced, hence increasing the drain current.. 

Next, consider the three terminaX FET to be operated: in-, 
the saturation region with a constant gate-source eraf.. There:-- 
fore, Ipg = The differential equaition. fair the drain-, 

current is given [Ref. 8 ] by the following: relationship:; 



dIpg/dT = Olp3/3i|;)(di[»/dT) + (3Ij^3/3a) (do/dT) 



= g^ (di|;/dT) + (Ijj 3 /o)(da/dT) 



’m 



4.46 



since 91 ^^ 3 / 91 !; = g^^ and 1^3 is proportional to a at a given 

\ 

temperature T. A reasonable approximation to the mobility- 
temperature function [Ref. 2] is 

y = p^(T/T^)-n 



o " O' 

where n depends on the impurity concentration. Therefore, after 



necessary manipulation the following equation is obtained: 

(l/o) (do/dT) = -n/T 4.47 

If AT << T then dip/dT can be approximated by ±2.2 mV/°C 
(depending upon whether the transistor is p-channel or 
n-channel, respectively) at low source curren.t Ieve:L.. For many 
purposes n may be approximated to 2 (actually it ranges from 1.7 to 
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2.7). Taking the above appro 2 ±inaEtr±(ms iirnGo: aci:iniirtt ttiee 
Eq. 4.48 can be obtained far ducEiini caiarceart: GhaargEe w.iitf. 
respect to temperature. 

dIpg/dT = 2.2(mV/"C) - n L^/T 4-. 45 

It is easily seen from Eq. 4.4^S tJiait: tdie: cimdiitlionr. fbxr zero: 
temperature coefficient is appra^cdnrateLy gxv'.err by 

W§m " 2.2 (T/h) mV .. 44 45? 

When the change in temperature is very smalX compared: to: 
the given temperature, T Eq.. 4^..4'S carr bo:- expressed 

in the following form: 

AIpg/AT = 2.2(mV/"C) g^ - n. L^^/T 

-then 

AIj ^3 = 2.2 g^ (AT) - (n/T) Ij^s(AT) 4-. 5D 

4 

Therefore, one may represent Eq. 4.50 using two current 
generators with current directions opposing to each other. 
Applying the Eq. 4.50 to Fig. 4-17 the FET computer model 
is obtained for a given temperature T and the temperature 
change AT << T, as shown in Fig. 4-25. If AT is^ la.rge, one 
has to make iterative computations by taking proper AT 
increments . 

The necessary model parameters may be determined as 
explained in the previous section. The currents and I 2 
caused by the temperature can be expressed in terms- of known 
parameters as 

Ij = 2.2 (AT) 4;51'. 
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l 2 = (n/T) Ijj 3 (AT) 4.52 

To be dimensionally correct g should be in millimhos, 

c>m » 

temperature T and (AT) in degrees Kelvin. 




Figure 4-25. Lumped JFET Computer Model with Temperature 

Effect . 

The enhancement- type MOSFET has a characteristic 
threshold voltage, V^, instead of the pinch-off voltage, Vp, 
but the equations are of the same general form: for 0 < v^^g 
- < Vpg (pinch-off region) i^g = K(Vgg - V,j,) . Here K is 

a constant and is a function of channel geometry, doping, and 
effective carrier mobility. In this section an analysis 
is given for the temperature dependence of Vgg when the 
current 'is maintained constant and saturation conditions 
are present. It is assumed that interface charge Qgg is 
independent of temperature and that the effects of a finite 
drain resistance can be neglected. Further, it is assumed 
that the only temperature dependent quantities are Vgg (DC 
gate-source voltage) , (Fermi potential) and (effective 
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surface mobility) . The temperature coefficient as a 
function of is given [Ref. 8] by 

dVcs/dt = (ip/gJ(-l/p^)(dy^/dT) 

4.53 

- (dc[)p/dT){[-2(V(.3 - V^)/Vjjg] + 1} 

Since the first and the second terms are both positive, 
it is evident that there may exist a realizable zero tempera- 
ture coefficient point. Supposing that the mobility is pro- 
portional to T and letting n = 2, one can obtain (by 

taking certain known factors and coefficients into account) 
the following expression for temperature coefficient [Ref. 8] . 

dvgg/dt = (3.0/T)(V(,3 - V^) 

- [(0.605 - (f)p)/T]{[2(V(,3 - V^)/Vp3l -1} . 

As pointed out by Giralt et al . [Ref. 23] , ‘a good approxima- 
tion to the behavior is given by 

dIvggl/dT = h - 2.5|v^3 - V^|mV/°C , 4.55 

where V.p is the measured threshold voltage obtained from the 
straight line extrapolation of an plot on the v^g axis, 

and h is a constant whose value is somewhat dependent on the 
subsrate doping but normally lies between 2 and 6. 

If AT << T then the Eq. 4.55 can be written as in 
Eq. 4.56. 

Ajvggl = h(AT) - 2.5|vgg - V^|(AT)mV . 4.56 

On the other hand when temperature changes as much as AT 
Eq. 4.43 becomes 
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^DST = ^ ' ^GS^ ^ ^DS 



" ’^GS "■ ^ml^^Gsl ‘ ^m ^GS ^DS 
" ^DS ^ml^^Gsl • 

This change Av^g can be obtained either from Eq. 4.54 or 
Eq. 4 .56 depending upon the given (or known) information. 
The term (gjj^ I I ^ defined as temperature resultant 

current generator than temperature dependent computer model 
for MOSFET obtained as shown in Fig. 4-26. 




Figure 4-26. Lumped • MOSFET Computer Model with Temperature Effect 
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V. coMCHJsrotfS 



A. GENERAL 

It has been shown that temper atujre \rar±at±(ins? hav.ee 
various effects on different electramic components;-.. Resrse- 
tance of a resistor increases or decreases w.iithe tera.pecr.atiir:e. 
depending on the material used. Inductance of an. iindiictnor: 
may or may not change with temperature depending; on its physii-- 
cal dimensions and the material used. Capacitance of: a 
capacitor changes slightly with temperature.. However:-,, this.- 
small change can be extremely important in. scone siitnationsi 
(e,g., as in resonance circuits). Forward diode currents- an cL 
diffusion capacitance increase with increasing temperature. 

4 

The effect on a transistor is quite complex, and some simpli- 
fying assumptions are listed in Chapter IV. The junction 
forward current and diffusion capacitance increased with 
increasing temperature. Reverse saturation current also 
increased. Two types of field-effect transistors were 
examined separately. It was found that with the proper 
operating point, the temperature effect could be eliminated. 

If this operating point cannot be attained, the drain current 
changes. Therefore, the effects discussed above should be 
applied to the computer models when one attempts to analyze 
a circuit which is composed of electronic components. 

The models proposed here allow simulation oi^" thu behavior: 
of several circuit components. For a reasonably accurate 
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simulation, the parajiietieirs ireaquiiiireidl finr thee nroaiEdss wearee few 
easily determined. Tiie atGrCuad: cxiixmiit; im wHliclr; thee conrpDnent: 
is to be used is a major faEdror- rir. thee dro±ae£ off ae mo dell andf 
the use of these models mus^t be done with. an. underst-airding of ' 
the assumed operating regioirs mrdL ch',amct3eihfsctdjt55 off thee 
components themselTes. The paraineteai^ off thee mocdedl nrastt also: 
be chosen to suit the particular app:Tlcaf±on , , berausee earenr, 
the best models can only approecimate the behavlorr off thee 
actual devices. 

Emphasis was placed on enrp^lrxcal defe.Tmihafii)n', off parameters ; 
however, these parameters^ may also: be deteoeminedf from: thee 
material characteristics and structural geometry off the: devices: 
being modelled. 

These models are useful in analysis and design of’ elec- 
tronic circuits. Also, since the model parameters may be 
calculated from information on the construction of the device, 
it should be possible to predict the behavior of a type of 
component before it is fabricated. 

The models are general so that they are compatible, with 
most of the electronic circuit analysis programs in use today. 

B. SUGGESTIONS FOR FURTHER STUDY 

1 . Modification of Models with High Frequency Effects 
High-frequency models for the components are more 
complex than low frequency models,, since additional, reactive 
elements must be included.. As= a suggesflon- for- further study,, 
one might consider the effects of high frequency effects 
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upon these components mideir mnmralL aqMa::atlirgr ttmpsTHcttiTESf 
and temperature as a Trariatrle.. 

2. FET in Triode Reg-ion 

As far as the FET is concerned,, the. developments 
were done for the saturation: region:, operation:.. Thee other: 



very important application is the- triadje-region-: ope-rationr. , 



In triode region, 
as 



Ca < « « 






- Vp-J' 



the- g _|)25 e-xpressedc 






5ii: 



At an arbitrary point oiitputr characteristic: 

curves, the pointwise linear approximation: to: thee output: curve 



is the straight line (Vj^g- - V^g-} 



where; gps- 



is the slope at (Vj^g, Inc) g'iven with Eq.. 5.M1.. Then' it: 



DS 



should be simple to develop models for triode region for 
various FET’s by taking the temperature effects into account. 

3 . Some Other Device Models 

One might also develop computer models for additional 
devices depending upon the methods stated in this- paper' with 
temperature or high frequency effects in mind. 

(a) Zener/avalanche diode madel: Since the avalanche 
diode is nearly identical to an ordinary diode for 
operating voltages greater than Vg (breakdown voltage) , 
it should be possible to modify the Fig. 4-9, diode 
computer model to include the avalanche phenomenon. 

(b) Tunnel diode model:: The unique properties: of: tunnel . 

diodes which are most often utilized are its negative 
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resistance and its very rapid switching from Vp 
(peak voltage) to a point in the forward conduction 
region. Most tunnel diode circuits are designed around the 
peak point and valley point. Any tunnel diode model should 
accurately approximate these points as well as modelling 
other important properties. It is relatively easy to 
measure the peak point, valley point and the maximum 
slope in the negative conductance region. There is 
very little resistance to current flow when reverse bias 
is applied to the junctions; i.e., with reverse bias, 
the tunnel diode is nearly a short circuit. 

(c) Unijunction transistor: On the emitter characteristics 

the peak point is defined 

< ^B2 




by the ex;f)ression 

V = V + aVtjn where v 
p BB 

is defined by the 1^2 “ ^ 
curve. The quantity a, 
called the intrinsic 

standoff ratio. When Ip becomes greater than 1 ^ (peak point 
current) , the process of conductivity modulation occurs in 
the semiconductor material between emitter and base- 1 , 
causing the negative resistance characteristic. As Ip is 
increased well beyond 1 ^ there is also a slight conductivity 
modulation of the - B2 region resulting in an increased 
Ip2‘ A model which will adequately describe the static 
behavior of a unijunction transistor can be formulated. 
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(3) Silicon controlled rectifiers:: Sinrcee SFfllii 

rectifier (SCR) is a facir-Iayery tJiree— junirt±an: dia/'.irxe w.itti-', 
only three leads, it is difficult to descnilie the irrtexnall 
behavior. The two-transistor analogy pTOviides a^ g.aojd: 
intuitive feel for the aperation,. frut iit rss daifSdiciiljt: tu: 
obtain the transistor parameters since bnse— h and: CDiIljsE±orr2 
are not accessible. It should be possihle toe repnessent: 
an SCR by empirically relating the internal behaviox- to: 
terminal characteristics even the physical processes' are 
somewhat more complex than for some of the devices^ dls^- 
cussed previously. 
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